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ABSTRACT 

We report on the sensitivity of SPIRE photometers on the Herschel Space Observatory Specifically, we measure the confusion noise from ob- 
servations taken during the Science Demonstration Phase of the Herschel Multi-tiered Extragalactic Survey Confusion noise is defined to be the 
spatial variation of the sky intensity in the limit of infinite integration time, and is found to be consistent among the diff'erent fields in our survey 
at the level of 5.8, 6.3 and 6.8 mJy/beam at 250, 350 and 500 //m, respectively. These results, together with the measured instrument noise, may be 
used to estimate the integration time required for confusion limited maps, and provide a noise estimate for maps obtained by SPIRE. 



Key words, cosmology: observations — confusion limit — infrared: diff'used background 



1. Introduction 

The Spectral and Photometric Imaging Receiver (SPIRE) 
(IGriffin et al. 20101) onboard the Herschel Space Observatory 
dPilbratt et al. 2010t has opened a new window on the Universe 
at far-infrared (FIR) wavelengths. The sensitivity of SPIRE de- 
tectors combined with HerscheFs 3.5 meter aperture allow as- 
tronomers to observe the EIR sky with unprecedented efficiency. 
With 18.1", 24.9", 36.6" (EWHM) beams at 250, 350 and 500 
microns, respectively, we expect SPIRE maps to be dominated 
by confused sources. It is therefore useful to determine the key 
characteristics, both of the instrument and the sky, that would al- 
low observers to optimize their observing plans, and/or to make 
sky surveys that probe as deep and as wide as possible for al- 
located observing time. In this letter we report measurements 
in SPIRE maps of the following: i.) instrument noise, i.e. noise 
from the detectors, readout electronics and photon noise from 
the telescope, ii.) confusion noise, that is, the variance in the sky 
map due to the presence of unresolved sources, and iii.) cosmic 
variance, arising from underlying large-scale fluctuations in the 
galaxy number density. 

The confusio n noise due to sources below a given flux cutofl', 
Xc, is derived in Condon 19741 to be the second moment of the 
measured flux distribution: 
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* Herschel is an ESA space observatory with science instruments 
provided by by European-led Principal Investigator consortia and with 
important participation from NASA. 



where x is the measured flux, x = S f(6, 0), S is the source flux 
convolved with the normalized beam response, f(6, cp), and dn 
is the difl'erential source distribution. In this analysis, the mea- 
sured fluxes have been binned into sky map pixels and CTconf cor- 
responds to the variance in a sky map in the limit of zero instru- 
ment noise, or equivalently, infinite integration time. This physi- 
cal definition allows us to derive a direct, simultaneous measure- 
ment of both the confusion noise, CFconf, and instrument noise, 
(Tinst^ which then allow us to define the SPIRE confusion limit. 

In Sect. 2, we describe the observations, followed by the 
analysis details and results in Sect. 3. In Sect. 4 we provide 
the number of "repeats" required to achieve confusion limited 
maps. A repeat, as defined by the ''Herschel Space Observatory 
Observation and Proposal Submission Tool" or HSPOT, contains 
two orthogonal scans of the entire field. At nominal scan speed, 
a single repeat averages to 16, 29 and 28 samples/pixel, or an 
integration time of 0.9, 1.6 and 1.5 seconds/pixel. In addition, 
we discuss the often-used (albeit, model-dependent) definition 
of confusion limit in terms of the source density or number of 
sources per 30 to 40 beams. 

2. Data Sets 

The Herschel Multi-tiered Extragalactic Survey (HerMESJj) 
Science Demonstratio n Phase (S DP) observations are detailed 
in Table 1 of Oliver e t al. 201 In this letter we have mea- 
sured the noise properties of SPIRE maps of three of the five 
fields, GOODS-N, Lockman-North and Lockman- SWIRE. The 
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Fig. 1. The GOODS-N field at 250, 350 and 500 yum (left to right). Only the central 16'xl6' of the 30'x30' field is shown above, 
with the green circles indicating the beam FWHM in each band. As expected, many sources can be seen at 250 yum that are not 
resolved at longer wavelengths due to the larger beams as many of the 250 jim sources would be resolved with a smaller beam. 
These images are a good illustration of confusion. 



GOODS-N field is 30'x30' and was covered by 30 map repeats. 
The Lockman-North field is 35'x35' covered by 7 repeats. These 
fields were observed in nominal scan mode with the space- 
craft scanning at 30"/s with an angle of ± 42.4 degrees with 
respect to the spacecraft y-axis (see SPIRE Observers Manual 
(2010)). Lockman-SWIRE is 218'x218' wide and was observed 
in SPIRE fast scan with a 60"/s scan speed, and was covered 
by 2 repeats. Standard SPIRE pipeline maps (see Oliver et al., 
in prep.) were used. These maps are calibrated in mJy/beam 
and have pixel dimensions of 6", 10" and 14" at 250, 350 and 
500 yum, respectively. The two remaining HerMES SDP fields, 
ELS and Abell 2218 were not used in this analysis. The ELS 
was contaminated by Galactic cirrus which is clearly seen in the 
measured pixel variance. Abell 2218 is a small cluster field with 
lensing eff'ects. Both of these fields are strongly biased to higher 
confusion noise. 

The central 16'xl6' of the GOODS-N field are shown in 
Eig. [T] at 250, 350 and 500 yum, clearly demonstrates the eff'ect 
of confusion in SPIRE maps. We note that the confusion noise 
according to eq. (1) is defined with some cutoff' ffux, Xc, which 
is set by the brightest sources in these maps. We will discuss 
the eff'ect of these cutoff's to the determination of the confusion 
noise. 



3. Analysis 

Eundamentally, spatial ffuctuations in a map arise due to two 
sources: instrument noise and the presence of sources on the sky. 
The instrument component, crinst, will be reduced with integra- 
tion time, while ffuctuations due to the convolution of the sky 
with the instrument beam, (Tco„/, will remain. Since the sources 
of noise are uncorrelated, and if low-frequency correlated noise 
has been properly accounted for in the low level data reduction, 
the total variance in the map is simply, cr^^^^/ = o-^onf "^ ^Instl^' 
where t is the integration time. Note that (Tconf has unit of surface 
brightness in mJy/beam, while crinst in mJy/beam ^[s. A straight 
line fft of the variance vs. inverse integration time has a slope 
that determines the instrument variance, cr^ , and a non-zero in- 



Table 1. Measured SPIRE Noise 



tercept that determines the variance of the sky intensity, cr^^„^ 
(IJarosik et al. 20031) . 

In our analysis we begin by selecting all map pixels with a 
total integration time between t - dt < t < t -\- dt, where dt is 



Band 



(^conf 

mJy/beam 



mJy/beam V^ 



250 
350 
500 



5.8 + 0.3 
6.3 + 0.4 
6.8 + 0.4 



8.5 + 0.4 
9.4 + 0.5 
13.3 + 0.7 



mJy/beam ^jrepeats 
Nomina^ Fasfl 



9.0 

7.5 
10.8 



fL 



12 

10.6 

15.3 



^ 30"/s scan speed 
^ 60"/s scan speed 



0.5 times the integration time per sample (i.e. 0.054 s/sample in 
nominal scan mode). We then measure the variance, cr^^^^/, of 
the pixel values in that subset. EigureS shows cr^^^^/(0 vs. r^ at 
250, 350 and 500 yum in the ffelds observed in nominal and fast 
scan mode. GOODS-N pixel sets are shown in black, Lockman- 
North in green and the shallow Lockman-SWIRE ffeld in red. 
The black line shows a simultaneous linear fft to all three ffelds 
and traces a component proportional to t~^ and a non-zero inter- 
cept. The data from the three included ffelds clearly show that 
the noise properties of the SPIRE instrument and the SPIRE ob- 
served sky are independent of the SPIRE scanning speed and sta- 
ble from shallow to deep ffelds. The consistency of the result in 
deep and shallow ffelds conffrms that this method of noise char- 
acterization measures the confusion noise in maps that, them- 
selves, are not confusion limited. In fact, the large area shallow 
ffeld (red points) in Eig. |2] show far less scatter about the fft than 
the fewer pixels in the deep ffelds since there are many more 
pixels, providing a larger statistical ensemble. 

The right side of Eig. [2] shows the same data in a diff'erent 
graphical representation. The square root of the variance in all 
three ffelds (now all in black) is plotted vs. the integration time 
along with the derived instrument noise (red line) and confusion 
noise (blue line). If the sky were dark, the noise would integrate 
down to zero following the red line, while for the real sky the 
noise can not get below the confusion ffoor. 

The measured noise is presented in Table 1 . The confusion 
noise values in Table 1 are equivalent to the \cr point source sen- 
sitivities, in mJy, in confusion limited maps. The two rightmost 
columns are included for observation planning and show the in- 
strument noise scaled from mJy/beam ^[s to mJy/beam Ajrepeats 
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Fig. 2. LEFT: Variance of SPIRE map pixels at 250, 350 and 500 jim (top to bottom) in 6", 10" and 14" pixels, respectively, in the 
GOODS-N (black), Lockman-North (green) and Lockman-SWIRE (red) fields. Lockman-SWIRE was only observed once in each 
scan orientation, and was observed in Fast Scan Mode (60"/s). Despite the shallow field depth and diff'erent scan speed, the linear 
trend with inverse integration time is clearly continuous from short to long integration times. The black line is a linear fit to all pixels 
in the three fields simultaneously revealing an instrument noise of 8.5 ± 0.4, 9.4 ± 0.5 and 13.3 ± 0.7 mJy/beam ^/s, and a confusion 
noise of 5.8 ± 0.3, 6.3 ± 0.4 and 6.8 ± 0.4 mJy/beam at 250, 350 and 500 yum. RIGHT: Pixel noise vs. integration time for all pixels 
in both fields. Overplotted are the derived instrument noise (red) and the confusion floor (blue) and the total noise (green). 



for nominal and fast scan mode, respectively, by averaging the 
total integration time per pixel in a map for a given number of 
map repeats. Fast scan mode is V2 higher since it has half the 
samples in a sky pixel for a given number of map repeats. 

We have made 10 simulated realizations of our survey fields 
by constructing timestreams with known instrument noise and 
injecting sources drawn from the number counts determined by 
BLAST (Patachon et al. 2009). These simulations are used to 
check for biases due to the non- gaussian nature of the pixel 
distribution or correlated noise and to provide better estimated 
uncertainties in our measurement. These simulated timestreams 
(Glenn et al., in prep.) are then processed through the mapping 
pipeline, and the above variance analysis is run on the simulated 
maps. The confusion and instrument noise measured in the sim- 
ulated maps agree with the analytically computed input noise in 
all three bands, indicating that any bias is insignificant. The con- 
sistency of our results in the deep and shallow fields further in- 
dicates an absence of measurement bias due to non-gaussianity. 
The scatter in the instrument and confusion noise determined 
in these simulations about the known inputs indicates that the 
statistical uncertainty on the instrument noise is 5% and on the 
confusion noise 6%, with no significant bias in the recovered 
values. 

The agreement among the three fields of varying depth and 
solid angle also indicate that an additional uncertainty due to 
cosmic variance is unnecessary. Cosmic variance, often negligi- 
ble in relatively large surveys, can become a significant source of 
uncertainties in deep pencil-beam surveys for high redshift stud- 
ies. As a check, we re-ran the variance analysis on sub-fields of 
the GOODS-N field to estimate the field size at which the noise 
estimates diverge from those made from the full dataset. That 
analysis shows that the confusion noise reported here is unaf- 
fected by cosmic variance in fields larger than 8'x8', lO'xlO' 
and 13'xl3' at 250, 350 and 500 yum. As our smallest map is 
30'x30', we include no uncertainty due to cosmic variance. 

Still, the accuracy of these results may depend on the data 
product or the quality of the map. For independent verification. 
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Fig. 3. Histogram of the 250 yum GOODS-N map, in black. The 
null map or jackknife removes the sky contribution and gives a 
map of the instrument noise, shown in the dashed histogram with 
the gaussian fit in red. The instrument noise is integrated down to 
1.7 mJy/beam, or 9.3 mJy/beam ^^[repeats. Similar analyses at 
350 and 500 yum give instrument noise of 1.4 and 1.9 mJy/beam 
or 7.7 and 10.4 mJy/beam ^/repeats, consistent with the values 
in Table 1 . 



further analyses have been performed using an independent map 
maker (SMAP iterative map maker, Levenson et al., in prep.). 
Two maps of the GOODS-N field were made, one using the 
first 15 repeats and one with the second 15 repeats. A null map 
or jackknife was then created by subtracting these two maps. 
Figure |3] shows the pixel histogram of the GOODS-N map along 
with the first/second half jackknife map at 250 yum. The jack- 
knife should remove sky signal and leave only instrument noise. 
Fitting a gaussian to the noise map histograms for GOODS-N 
maps with various numbers of repeats yields an instrument noise 
of 9.0 ± 0.1, 7.2 ± 0.1 and 10.4 ± 0.4 mJy/beam ^repeats at 
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Fig. 4. Variance plot for jackknife map of GOODS-N. As ex- 
pected, the null map gives a linear fit with the similar slope or 
instrument noise as the full map, and the intercept or confusion 
noise is consistent with zero since sky contribution has been re- 
moved. 
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Fig. 5. Confusion noise as a function of pixel flux cut at 
250/350/500 yum in blue/green/red. All pixels within a beam 
FWHM of a pixel brighter than the indicated flux cut were 
masked and the confusion noise was re-estimated. Vertical lines 
indicate ScTconf in each band. 



250, 350 and 500 yum. These values are consistent with those 
measured using the variance method described above (see Table 
1) and, as the variance method gives a precision measurement of 
the confusion noise as well, we present the results of that analy- 
sis as the fiducial SPIRE noise estimates. 

Figure |4] shows the equivalent of left column of Fig. [2] for 
the jackknife map of GOODS-N. As expected, the instrument 
noise remains in the jackknife map, and with the sources re- 
moved the intercept is now consistent with zero. Specifically, 
intercepts from jackknife map give noise estimates of 0.3 ± 0.3, 
-0.3 ± 0.3 and -0.4 mJy/beam ± 0.5 at 250, 350 and 500 jim. 
Instrument noise measured using the same method in the jack- 
knife map is 8.4 ± 0.9, 9.0+1.1 and 12.8 ± 1.8 mJy/beam V^ 
which agrees with the values determined from the full maps. 

4. Discussion 

We have made a precision measurement of the confusion and in- 
strument noise in the HerMES SDP SPIRE scan maps. The var- 
ious checks presented, including simulated datasets with known 
noise properties and a repeat of the analysis on null map or jack- 
knife, demonstrate the robustness of the analysis and give good 
consistency in the values of the confusion and instrument noise. 
Our simulations indicate that the statistical error of the confusion 
noise is about 5%. We note that the systematic error is largely 
dependent on the calibration and map-making process, currently 
estimated to be ~ 15% fGriflS n et al. 2010b . Table 1 can be used 
as a guide for achieving confusion limited maps with SPIRE. In 
particular, to make maps in which the instrument noise is compa- 
rable to the confusion, it will take a minimum of 3 map repeats in 
nominal scan mode or 5 map repeats in fast scan mode. Longer 
integration time will result in maps in which confusion is the 
dominant source of uncertainty in measurements of source flux 
and position. 

The noise values reported in Table 1 make no significant 
source cut and, accounting for the possibility that even bright 
sources are confused, measure the variance in HerMES maps 
due to all sources up to lOcFconf or 80 mJy (more than 99% 
of data). It may be interesting to systematically remove bright 
sources and study how confusion noise might change accord- 



ingly. We have determined the confusion noise using the same 
method presented in Sect. 3, after removing pixels within a beam 
FWHM of any pixel brighter than a given flux cut. (Since the 
maps are calibrated such that the source flux is given by the flux 
in the pixel at the peak of the PSF, i.e., mJy/beam, this method 
is equivalent to removal of sources. Removing pixels within a 
FWHM is conservative and ensures that any extended emission 
is completely removed.) The results of this analysis are shown in 
Fig. l5] The vertical lines indicate ScTconf as measured in the full 
map /The residual noise estimates are 3.8, 4.6 and 5.2 mJy with 
a ?>o-conf cut and 4.8, 5.5 and 6.1 mJy with a ScTconf cut. 

The confusion limit is often quoted in terms of the flux, 
at which the source density in a map reaches 30 to 50 beams 
per source. This quantity requires a measurement or a model of 
the density of sources, and often is derived assuming a power 
law source distribution with a 3 -Act con f JFranceschini 19821) . 
BLAST (Devlin et al. 2009) predicts confusion limits corre- 
sponding to 40 beams per source for SPIRE of 22, 22 and 18 
mJy at 250, 350 and 500 yum, respectively. The best source 
counts available at the SPIRE wavelengths are probably the 
HerMES counts themselves, as reported by Oliver et al. (2010). 
The HerMES counts reach 1 source per 40 beams at 19.1 ± 0.6, 
17.7 ± 0.6, and 15.1 ± 1.8 mJy at 250, 350 and 500 yum, corre- 
sponding to 3.29, 2.81, and 2.60(Tco«/, and in rough agreement 
with blast's prediction. 

Acknowledgements. SPIRE has been developed by a consortium of institutes led 
by Cardiff Univ. (UK) and including Univ. Lethbridge (Canada); NAOC (China); 
CEA, LAM (France); IFSI, Univ. Padua (Italy); lAC (Spain); Stockholm 
Observatory (Sweden); Imperial College London, RAL, UCL-MSSL, UKATC, 
Univ. Sussex (UK); Caltech, JPL, NHSC, Univ. Colorado (USA). This develop- 
ment has been supported by national funding agencies: CSA (Canada); NAOC 
(China); CEA, CNES, CNRS (France); ASI (Italy); MCINN (Spain); SNSB 
(Sweden); STFC (UK); and NASA (USA). 



References 

Condon, J.J., ApJ, 188, 279 (1974) 
Devlin, M., et al.. Nature 458, 737-739 (2009) 
Franceschini, A., AP&SS, 86, 3F (1982) 
Griffin, M., et al., this issue (2010) 
Jarosik, N., et al. , ApJ, 148, 29-37 (2003) 



H. T Nguyen et al.: HerMES: The SPIRE confusion limit 



Oliver, S., et al., this issue (2010) 
Patanchon, G., et al. 2009, ApJ, 707, 1750 
Pilbratt, G. L., et al., this issue ( 2010) 

SPIRE Observers Manual, HERSCHEL-HSC-DOC-0789, version TBD, Date, 
Website, ESA 



Jet Propulsion Laboratory, 4800 Oak Grove Drive, Pasadena, CA 

91109, USA 

e-mail: htnguyen@jpl . nasa . gov 

California Institute of Technology, 1200 E. California Blvd., 

Pasadena, CA 91 125, USA 

Infrared Processing and Analysis Center, MS 100-22, California 

Institute of Technology, JPL, Pasadena, CA 91125, USA 

Dept. of Physics & Astronomy, University of California, Irvine, CA 

92697, USA 

Institute for Astronomy, University of Edinburgh, Royal 

Observatory, Blackford Hill, Edinburgh EH9 3HJ, UK 

Laboratoire AIM-Paris-Saclay, CEA/DSM/Irfu - CNRS - Universite 

Paris Diderot, CE-Saclay, pt courrier 131, F-91191 Gif-sur-Yvette, 

France 

Astrophysics Group, Imperial College London, Blackett Laboratory, 

Prince Consort Road, London SW7 2AZ, UK 

Laboratoire d'Astrophysique de Marseille, OAMP, Universite Aix- 

marseille, CNRS, 38 rue Frederic Joliot-Curie, 13388 Marseille 

cedex 13, France 

Institute de Astrofisica de Canarias, C/ Via Lactea s/n, E-38200 La 

Laguna, Spain 

Department of Physics & Astronomy, University of British 

Columbia, 6224 Agricultural Road, Vancouver, BC V6T IZl, 

Canada 

Dept. of Astrophysical and Planetary Sciences, CASA 389-UCB, 

University of Colorado, Boulder, CO 80309, USA 

Herschel Science Centre, European Space Astronomy Centre, 

Villanueva de la Canada, 28691 Madrid, Spain 

Observational Cosmology Lab, Code 665, NASA Goddard Space 

Flight Center, Greenbelt, MD 20771, USA 

Cardiff School of Physics and Astronomy, Cardiff University, 

Queens Buildings, The Parade, Cardiff CF24 3AA, UK 

Dipartimento di Astronomia, Universita di Padova, vicolo 

Osservatorio, 3, 35122 Padova, Italy 

ESO, Karl-Schwarzschild-Str. 2, 85748 Garching bei Miinchen, 

Germany 

UK Astronomy Technology Centre, Royal Observatory, Blackford 

Hill, Edinburgh EH9 3HJ, UK 

Institut d'Astrophysique Spatiale (IAS), batiment 121, Universite 

Paris-Sud 11 and CNRS (UMR 8617), 91405 Orsay, France 

School of Physics and Astronomy, The University of Manchester, 

Alan Turing Building, Oxford Road, Manchester M13 9PL, UK 

Astronomy Centre, Dept. of Physics & Astronomy, University of 

Sussex, Brighton BNl 9QH, UK 

Institut d'Astrophysique de Paris, UMR 7095, CNRS, UPMC Univ. 

Paris 06, 98bis boulevard Arago, F-75014 Paris, France 

Mullard Space Science Laboratory, University College London, 

Holmbury St. Mary, Dorking, Surrey RH5 6NT, UK 

Space Science & Technology Department, Rutherford Appleton 

Laboratory, Chilton, Didcot, Oxfordshire 0X11 OQX, UK 

Institute for Space Imaging Science, University of Lethbridge, 

Lethbridge, Alberta, TIK 3M4, Canada 

Astrophysics, Oxford University, Keble Road, Oxford 0X1 3RH, 

UK 

Centre for Astrophysics Research, University of Hertfordshire, 

College Lane, Hatfield, Hertfordshire ALIO 9AB, UK 



